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Abstract 

We derive self-reciprocity properties for a number of polyomino generating func- 
tions, including several families of column-convex polygons, three-choice polygons and 
staircase polygons with a staircase hole. In so doing, we establish a connection be- 
tween the reciprocity results known to combinatorialists and the inversion relations 
used by physicists to solve models in statistical mechanics. For several classes of con- 
vex polygons, the inversion (reciprocity) relation, augmented by certain symmetry 
and analyticity properties, completely determines the anisotropic perimeter generating 
function. 
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1 Introduction 



Symmetries are among the most important guiding principles in all of physics and mathe- 
matics. It often happens that a problem may be solved by symmetry considerations alone, 
and even if not, understanding the symmetries of the solution can greatly reduce the amount 
of work needed to find it. We study here a symmetry of functions which is known as "self- 
reciprocity" to combinatorialists and which is referred to as "inversion relations" in lattice 
statistical mechanics. 

Our focus will be on polyomino enumeration problems which are of interest in both 
combinatorics and physics. We shall demonstrate that one can find examples of functional 
symmetry in the resulting generating functions. 

The inversion relation rose to prominence in statistical mechanics in the early 1980s as 
the most direct path to the solution of many integrable models |^3|, 0, [| and was soon 
realized to be commonplace in both solved and unsolved models [||, |3], Let G(x) be 
a thermodynamic quantity which depends on a collection of parameters, x. An inversion 
relation is a functional equation 

G{x) ± x a G{(f){x)) = ifj{x) (1) 

where <f> and if) are known functions of x. Typically (f> involves taking reciprocals of one or 
more components of x. The inversion relation tightly constrains the function G. For some 
two-dimensional models a pair of additional conditions holds: that G is symmetric under 
exchange of horizontal and vertical, and that G is an analytic function of its arguments. 
Very often, the three constraints taken together uniquely determine the function G. 

In 1974 Stanley presented a general framework for reciprocity results. He established sev- 
eral powerful general conditions under which a generating function will be self-reciprocal 
The language and notation of Stanley p0|, ^2] will be used throughout this paper. 



Definition 1.1. Let H(yi, . . . , y n ) be a rational function in the variables yi, with coefficients 
in R. We say that H is self-reciprocal if there exists an n-tuple of integers (fix, . . . , (3 n ) such 
that 

H(l/y h ... ,l/y n ) = ±y? 1 ...y^H(y 1 ,... ,y n ). (2) 



In what follows, we write y& = yf 1 . . .y@ n and 1/y = (1/yi, ■ ■ ■ , l/y n )- Thus eqn. (H) may 
be concisely expressed as 

H(l/y) = ±y< 3 H(y). 
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Note that a rational function is self-reciprocal if and only if both its numerator and denom- 
inator are so, and that the self-reciprocity of a polynomial amounts to a certain symmetry 
in its coefficients. Some explicit examples are given in Subsection |3.2| . 

Let us now demonstrate the relationship between self-reciprocity and inversion relations. 
Consider the multivariable generating function 

G(x, y) = J2 C(m, n)xT • . • x?y? . . . y n k k 

= J2C{m,n)x m y n (3) 

where m = (mi, . . . , rrij), x = (x lf . . . , Xj) and similarly for n and y. The summation is over 
(j + fc)-tuples of nonnegative integers representing the objects being enumerated. Performing 
the summation over n, we reexpress eqn. in terms of partial generating functions, H m (y), 

G(x,y) = Y J H m (y)x m . (4) 

m 

Now suppose that the partial generating functions are self-reciprocal, 

H m (l/y) = ±e m y^H m (y), (5) 

where e is a j-tuple of elements in the set { — 1, 1} which characterizes the dependence of the 
sign on m, and where /3(m) depends linearly on m: 

f3(m) = Am + a. (6) 

Here, A = (a^i)^, is a k x j matrix of integers and cx is a fc-tuple of integers. We can then 
write 

G(x, y) T y a G(exy A , l/y) = (7) 

where exy~ k is the j-tuple whose i th entry is e^Xj n^=i Vt 0< '*- This is clearly a special case 
of the inversion relation ([!]). A few comments are in order: 

• The right hand side of (0) is zero, but in the more general situation some of the partial 
generating functions, H m (y) will fail to be self-reciprocal for certain choices of m. If 
we are fortunate, this will be a small, finite or otherwise controllable set of cases, and 
we will be able to compute the correction term we need to add to the right hand side 
explicitly. For many examples in statistical mechanics, this correction term depends 
on x but not on y. 
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In all of the cases we shall see below, the denominators of our rational functions will 
be a product of terms (1 — y aj ), which are self- reciprocal. Stanley has proved that this 
denominator form always holds for certain classes of problems (see Theorem 4.6.11 of 
ref. O). 



• It might be asked which of the concepts, inversion or self-reciprocity, is the more 
general. On one hand, in the derivation of (|7|) the dependence of the exponent f3 on 
m was assumed to be linear, which may not always hold, implying that reciprocity is 
more fundamental. On the other hand, the function occurring in ([I]) may in principle 
be more complicated than x — > exy A , y —>■ 1/y. In this case, the partial generating 
functions might not be self-reciprocal. An example is provided in Section by the 
Potts model, but in the polyomino examples considered in this paper, this situation 
does not arise. 

We now present a nonexhaustive list of recipes for finding and proving reciprocity results 
and inversion relations. 

1. If the generating function (or thermodynamic quantity) is known in closed form, an in- 
version relation can be demonstrated directly. As an example, we treat the anisotropic 
perimeter generating function for directed convex polygons in this manner in Section || 



For statistical mechanics models which admit a formulation in terms of a family of 
commuting transfer matrices, a transformation of parameters can often be found which 
inverts the transfer matrix. The commutativity property then allows the inversion 
relation to be derived. We review this in detail in Section 0, with the two-dimensional, 
zero- field Ising model as primary example. 

In non-integrable models, the transfer matrix will still be invertible and may suggest a 
possible inversion relation, but the required analyticity property is lacking. Neverthe- 
less, the suggested inversion relation can often be verified by inspection of the partial 
generating functions up to some finite order in the low-temperature expansion (Q). The 
q > 2 Potts model inversion relation discussed in Section p] was derived this way in 



ref. [13]. Some of the new results reported in the present paper were initially discovered 



by this method before being rederived by one of the other methods. 
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4. The "Temperley methodology" can be used to obtain very general reciprocity results 
for many classes of column-convex polygons. The first step is to derive a functional 
equation for the generating function which can be interpreted as the gluing of an 
additional column onto the graph. Step two is to show by induction that appending 
an additional column preserves self-reciprocity. This is detailed in Section [|. 

5. If the problem can be posed as a system of linear diophantine equations, whose solutions 
are subject to certain types of constraints, we may apply self-reciprocity theorems due 
to Stanley |p0[ . We have so far succeeded in applying this method only to families 
of directed polyominoes (Section [5]), but it enables us to treat problems which are 
impossible, or at least extremely cumbersome, by the method of functional equations. 



6. For combinatorial objects with a rational generating function of denominator fX/(l ~ 
t/"-?), one can try to explain self- reciprocity - i.e., the symmetry of the numerator - by 
interpreting the numerator combinatorially. This has been done by Fedou for a family 
of objects related to (but distinct from) staircase polygons JTOj . 

In Section |2] we review the motivation for looking at inversion relations in statistical 
mechanics and describe the methods used to obtain them. This will be useful for making 
comparisons with the results obtained later, and for suggesting applications and general- 
izations of the inversion relations. In Section |3|, we present examples of reciprocity results 
and inversion relations for polyominoes, and summarize our main new results. The technical 
heart of the paper consists of Section f|on the Temperley methodology, and Section [5] on the 
application of Stanley's results to polyominoes. 



2 Inversion relations in statistical mechanics 

The first use of the inversion relation in statistical mechanics was the solution by Stroganov 
of certain two dimensional vertex models on the square lattice p3 . Generalizations of 



Stroganov's models were later solved by the same means by Schultz |y|. Shortly after 
Stroganov, Baxter used a similar method to solve the hard hexagon model and rec- 
ognized its broad applicability, giving the eight- vertex and Ising models as examples [@. 
Subsequently, a number of authors pointed out that many known solutions to problems 
in two-dimensional statistical mechanics can be derived easily using the inversion relation 



method. Among these were Shankar [18], Baxter || and Pokrovsky and Bashilov [17 |. 
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It is noteworthy that inversion relations hold also for models that have not been solved. 
Prominent among such models are the two-dimensional Ising model in a magnetic field 
whose inversion relation was found by Baxter j3J, and the three-dimensional Ising model and 
noncritical g-state Potts model, both of whose inversion relations were found by Jaekel and 
Maillard [12, O. What generally distinguishes solved and unsolved models is the growth 



rate in the number of poles arising in the partial generating functions in the expansion 
as a function of order. Roughly speaking, a more complicated pole structure implies that 
the number of parameters needed to specify a given partial generating function is greater, 
and makes it less likely that an inversion relation can completely determine all of them. 
Nevertheless, inversion relations are still invaluable in the study of such problems, not least 
because they provide an independent check on series data. 

The fundamental problem of statistical mechanics is to calculate the partition function. 
Here we consider vertex models defined on a square lattice with each bond colored with one 
of r possible colors. Each lattice site makes a contribution to the energy of the system which 
depends on the colors of the adjacent bonds. This defines an r 4 - vertex model if all possible 
colorings are permitted. 

Stroganov computed the partition function per site in the thermodynamic limit of several 
16- and 81-vertex models. Consider first a finite lattice (on the torus) of iV rows and M 
columns. The partition function can be expressed in terms of the transfer matrix Tm as: 

Zm,n = Tr (Tm)^ 

(see ||, |22j). Here, T m is the r M x r matrix whose i, jth entry is the contribution to Zm,n of 
a single row of M sites connected to the row below by a set of vertical bonds in configuration 
i and to the row above by a set of vertical bonds in configuration j. It depends on the 
temperature, T, and on r 4 parameters specifying the vertex energies. In the thermodynamic 
limit, the partition function per site is given by 

k = Mm (Z M , N ) 1/MN = lim (A M ) 1/M (9) 

where Xm is the largest eigenvalue of T M , assumed to be nondegenerate. 

For simplicity let us consider a family of models whose vertex energies are functions 
of a single parameter, b. The models solved by Stroganov are integrable by virtue of the 
commutativity of the transfer matrices at different values of this parameter. This implies 
that the transfer matrix eigenvectors are common to all members of the family, and that 
the b dependence is only in the eigenvalues. For this reason b is often called the spectral 
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parameter. The key observation is that the inverse of the transfer matrix in these models is 
itself a member of the commuting family, up to a scale factor 

[XwP)]- 1 = T M (0(6)). (10) 

Acting on the eigenvector corresponding to Am (b) with both sides of eqn. (|T0|) yields the 
functional equation 

K(b)n( ( f ) (b))=m- (11) 

It is the commutativity of the transfer matrices for all values of b that allows the analytical 
continuation of the function k from b to <f)(b). With knowledge of the functions ip(b) and 
4>{b) and using the analyticity of n(b), Stroganov finds a unique solution, thereby reproducing 
Baxter's results for the symmetric eight vertex and homogeneous ferroelectric models, and 
obtaining the result for a certain 81- vertex model . 

As an illustrative example, we review here the derivation by Baxter of Onsager's 
expression for the partition function of the two-dimensional zero-field Ising model \TE\. Let 
the square lattice be drawn at 45° to the horizontal and let the couplings between nearest 
neighbors along the two lattice directions be J and J'. Define low temperature variables 

x = e ~ 2K } y = e~ 2K ' with K = J/k B T, K' = J'/k B T. (12) 

Transfer matrices for different choices of parameters will commute provided they have the 
same value of k = (sinh2i^ sinh 2K')~ 1 . The transformation which inverts the transfer 
matrix is 

K^K + Y , K'^-K', (13) 

which does not modify the value of k. Define the reduced partition function per site by 

A(x, y) = exp(-K - K')k(K, K'). (14) 

Then A(x,y) obeys the inversion relation 

A(x, y)A(-x, l/y) = 1 - x 2 . (15) 

Note that \ogA(x,y) — |log(l — x 2 ) has an inversion relation of precisely the form (|7p. 
By the symmetry of the model, we have 

A{x,y)=A{y,x). (16) 
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Inspection of the low temperature expansion leads us to conjecture the form 

m>l V y ' 

That the coefficient of x 2m is a rational function of y 2 is apparent from the nature of the low 
temperature expansion, but that the denominator has such a simple form is not expected 
on general grounds. Presumably it is a consequence of the condition of commuting transfer 
matrices. Here we take it as a hypothesis. Then Baxter has shown that the inversion 
relation fll5|) , symmetry (jlBf) and the denominator form ([i~7|) determine A(x,y) completely. 



We present his argument in Section where we use it in the context of polygon enumeration. 

Up till now we have been assuming integrability and in particular we have relied on the 
property that the transfer matrix and its inverse are both members of some one-parameter 
commuting family. What about models for which this property doesn't hold? Since analyt- 
icity of K,(b) breaks down, the step (pTTJ ) in the above derivation is no longer valid. However, 
it is still possible to obtain an inversion relation by direct analysis of the low-temperature 
expansion of the partition function to some finite order. As an example, it was shown in 



ref. that the logarithm of the reduced partition function per site, G(x,y) = lnA(x, y), 

of the g-state Potts model satisfies the inversion relation 

, „/ x 1\ f (1 - x)(l + (q - l)x) 

G(x,y)+G[--— — ,- =ln' 



l + (q-2)x , yj V l + (q-2)x 
When q = 2 this reduces to the Ising model inversion relation (|I5D . The inversion relations 
we will be considering in the remainder of the paper are derived by analysis of the generating 
function (analogous to the low temperature expansion) and do not depend on the models 
being integrable. 

An additional new feature is seen in this Potts model example. Neglecting for the moment 
the nonzero right-hand-side of (|18|), which can be eliminated by a suitable redefinition of 
G(x,y), we notice that when q > 2 there is no longer an order- by-order cancellation of the 
partial generating functions as defined in (|J), but rather cancellation of combinations of 
partial generating functions of different orders. However, we may convert to self-reciprocal 
form by defining 

^'''^( i-h-W ') (19) 

under which the inversion relation becomes 

G'{x, y) + G\-x, l/y) = In ^^^^ (20) 
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In the cases we will look at in this paper, the partial generating functions turn out to be 
self-reciprocal in the natural variables of the problem. We have not investigated the existence 
of inversion relations involving more complicated changes of variables. 

3 Polyomino enumeration and self-reciprocity 
3.1 Definitions 

The constructions we will consider are defined on the square lattice. All are defined only up 
to translation on the lattice. Starting at a lattice site and moving to one of the four nearest 
neighbors constitutes a step which we may identify with the edge connecting the sites. A 
connected sequence of steps is a path or walk. If no lattice site in the path occurs more than 
once, the path is self- avoiding. If a path returns to its starting site in the final step, and 
otherwise does not intersect itself, the result is a self-avoiding polygon. The number of steps 
taken is the perimeter of the polygon; the number of steps taken in the vertical direction is 
the vertical perimeter. The horizontal perimeter is defined similarly. The area is the number 
of cells of the lattice enclosed by the polygon. 

Enumerating self-avoiding polygons according to perimeter or area is an unsolved prob- 
lem. However, progress has been made in enumerating certain subclasses of self-avoiding 
polygons. Rectangles coincide with the rectangles of ordinary geometry whose vertices are 
lattice points and whose edges lie along lattice directions. A rectangle which contains a 
given polygon, i.e., all steps of the polygon lie inside or on the rectangle, is a bounding 
rectangle for that polygon. The smallest such rectangle is the minimal bounding rectangle. 
A polygon whose perimeter equals that of its minimal bounding rectangle is convex. If a 
convex polygon contains at least one of the corners of its minimal bounding rectangle (for 
concreteness say the south-west corner) then it is a directed convex polygon. If it contains 
also the north-east corner, it is a staircase polygon, so called because it is bounded above and 
below by two staircase-like or directed paths. On the other hand, if it contains two adjacent 
corners, say the southwest and southeast (northeast and southeast) then it is a stack polygon 
with horizontal (vertical) orientation. If it contains three corners, then it is a Ferrers graph. 
Representative examples of different classes of convex polygons are shown in Figure [I]. 

One way to obtain non-convex polygons is to relax the convexity condition along one 
direction only. A self-avoiding polygon is column-convex if the intersection of any vertical line 
with the polygon has at most two connected components. Row-conv ex polygons are similarly 
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(a) Ferrers graph 



(b) Stack polygon (horizon- (c) Stack polygon (vertical) 
tal) 



(d) Staircase polygon 



(e) Directed convex polygon 



(f) Convex polygon 



Figure 1: Classes of convex polygons 
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(a) Column-convex polygon 



(b) Bar-graph 



(c) Directed column-convex 
polygon 



Figure 2: Classes of column-convex polygons 



defined. The set of convex polygons is the intersection of the sets of row- and column- convex 
polygons. Subclasses of column-convex polygons include the bar-graphs which contain the 
bottom edge of the minimal bounding rectangle, and directed column- convex polygons whose 
bottom edge is a directed path. Some examples are shown in Figure 0. 

A second class of non-convex polygons is made up of four directed paths. A three-choice 
walk is a self-avoiding walk whose steps are taken in accordance with the three-choice rule 
which allows a step either to the left or the right or straight ahead after any vertical step, 
but forbids a right turn after any horizontal step. A polygon formed from such a walk is a 
three- choice polygon. When the walk returns to its starting point, we don't specify whether 
the next step, i.e., the first step, is a valid continuation of the walk. If it is, the result is 
a staircase polygon; if not, it is an imperfect staircase polygon (see Figure |3(a)|) . When we 
refer to three-choice polygons below, we include only the imperfect ones. 

A polyomino is a union of connected (sharing an edge) cells of the lattice. We shall 
consider one class of nonpolygon polyominoes — the staircase polygons with a staircase hole. 
The outer boundary and the hole are both staircase polygons and must not touch at any 
point. An example is shown in Figure 3(b) . 
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(a) Three-choice polygon (b) Staircase polygon with a (c) Self-avoiding polygon 

staircase hole 

Figure 3: Non-convex polyominoes 
3.2 Self-reciprocity in polyomino enumeration 

For each of the above classes of column-convex polygons, the anisotropic perimeter and area 
generating function, 

G(x, y, q) = E E C ^ n ' ^"W (21) 

m>l n>l a>l 

has been computed exactly (see ref. and references therein). Here C(m, n, a) is the number 
of polygons of the class with 2m horizontal bonds, 2n vertical bonds and area a. For the 
classes of convex polygons, the anisotropic perimeter generating function, G(x, y, 1) is an 
algebraic function of the fugacities, x and y, whereas the area generating function, G(l, 1, q) 
is a g-series. For classes of polygons that are only column-convex, both G(x, y, 1) and 
G{\, 1, q) p4| are algebraic, but G(x, y, q) involves g-series. A closed-form expression for the 
three-choice polygon anisotropic perimeter-area generating function is not yet known, but by 
means of a transfer matrix technique it can be evaluated in polynomial time ||. The isotropic 
perimeter generating function, G(x, x, 1) is known to have a logarithmic singularity 0, and 
is therefore not algebraic, but is known to be D-finite. The generating function for staircase 
polygons with a staircase hole is also not known in closed form. Its properties are expected 
to be similar in many respects to the generating function for three-choice polygons ||11|| . 

We shall be concerned with self-reciprocity properties of the generating functions H m (y, q) 
that count polygons of width m. We first give two examples. 

1. The area generating function for staircase polygons of width 4 is the following rational 
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function ||: 

g 4 (l + 2q + Aq 2 + 6g 3 + 7q 4 + 6g 5 + 4g 6 + 2q 7 + q 8 



(l- ? )2(l_ g 2 ) 2( 1 _ g3 )2 (l _ ?4) 

It satisfies 

i? 4 (l/g) = -H,{q), 

and is thus self-reciprocal. Observe that the numerator is not only symmetric (due to 
self- reciprocity) , but also unimodal. 

2. The (half-)vertical perimeter and area generating function for column-convex polygons 
of width 3 is the following rational function, which can be derived from the general 
formula of ref. flTfl: 



Hs(y, q) = Tl W1 m ^ =r ■ (yV + Ay 5 q 7 + 2y 5 q e + „V - yY 



f 

[i - m) 4 (i - m 2 ) 2 (i - yg 3 ) 

- Ay 3 q 5 - 6y 3 q A - Ay 3 q 3 - y 2 q 4 + y 2 q 2 + 2yq 2 + Ayq + 1). 



It satisfies 

H 3 (l/y,l/q) = ~H 3 (y,q) 

yq A 

and hence is self-reciprocal. Again, observe the symmetry of the coefficients in the 
numerator. 

We shall generalize these results to polygons of any width. Table [I] summarizes the self- 
reciprocity properties we have established. Most of them can be proved in various ways. 
One can for instance use a closed form expression of the generating function (Section |3.3| ), 
or a functional equation that defines it (Section f|); one can also encode the polygons by 
a sequence of numbers constrained by linear diophantine equations and apply Stanley's 
general results (Section ||). We shall see that the last two methods allow us to introduce 
many additional parameters and obtain self-reciprocity results that significantly generalize 
those of Table [p. 



3.3 Self-reciprocity via generating functions 

When a closed form expression for the generating function of some class of polygons is known, 
it seems natural to use it to demonstrate an inversion relation. Let us take the example of 
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Class 


Picture 


Self Reciprocity 


Inversion Relation 


Ferrers 

stack 
staircase 

directed convex 

convex 

bargraph 

dir. col.-conv. 

column-convex 




Li 
<\ 

£> 
ft 

h 


H m (l/y,l/q) = V ^ H m {y,q) 

H m (l/y,l/q) = -y 2m - 3 q m2 - 2m H m (y,q) 
H m (l/y, l/q) = -y m ~ l H m {y, q),m>2 

H m {l/y) = -y m ~ 2 H m (y) 

Not simple 

H m (l/y,l/q) = t^H m (y,q) 

H m (l/q) = -\H m {q) 

H m (l/y, l/q) = --L.H m (y, q) 


G(x,y)-y 2 G(-x/y,l/y) = 

G(x,y)+y 3 G(x/y 2 ,l/y) = 
G(x, y, q) + yG(x/y, 1/y, l/q) = -x 

G(x,y)+y 2 G(x/y,l/y) = 

G(x, y) + y 3 G(x/y, 1/y) = xy - x 3 y£^ 

G(x,y,q) -yG(-xq,l/y,l/q) = 

G(x,q) + qG(x,l/q)=0 

G(x, y, q) + yG(xq, 1/y, l/q) = 


three-choice 
SC with SC hole 






Not simple 
Not simple 


G(x, y, q) + y 2 G(x/ y, 1/y, 1/ q) = known 
G(x,y,q) + y 2 G(x/y, 1/y, l/q) = known 



Table 1: Summary of polyomino inversion relations 



the anisotropic perimeter generating function for directed convex polygons, which is known 
to be 0: 



G^y) = ^= (22) 

with A(x,y) = l — 2x — 2y — 2xy + x 2 + y 2 = (l — y) 2 [1 — x{2 + 2y — x)/(l — y) 2 }. Expanding 
this expression in x gives 



G(x, y ) = Y, H ^y) 



_ x 

m>l 

V , 3/(1 + 2/) 2 y(l + 4y + y 2 ) 3 y(l + 9y + 9y 2 + y 3 ) 4 5 
1-/ (1-1/) 3X (l-l/) 5 ^ (1-y) 7 x 

which suggests that the partial generating functions, H m {y) are self-reciprocal, and more 
precisely, that H m (l/y) = —y m ~ 2 H m (y). This is equivalent to the inversion relation 

G(x,y)+y 2 G(x/y,l/y) = J (23) 

which is easily checked from the closed form of the generating function. Note that an 
explicit expression for H m (y) is given in 0]. The inversion relations for convex polygons 
and directed column-convex polygons may also be obtained from the expression of their 
generating function. 

The partial generating functions for directed convex polygons, counted by the area, are 
not self-reciprocal: for instance, the generating function for width 3 is 

g 3 (l + 3q + 3q 2 + 2q 3 + g 4 )/(l - q) 2 (l - q 2 ) 2 (l - q 3 ). 

However, many other classes of column-convex polygons have an inversion relation for the 
full anisotropic perimeter and area generating function. Since these generating functions are 
also known in closed form they could be derived as above. However more can be shown, 
namely that there is a self-reciprocity for any parameter which is a linear function of the 
vertical heights in the graph. This very general result will be derived in Section |j. Likewise, 
the inversion relations for three-choice polygons and staircase polygons with a staircase hole, 
given in Table [I], are also special cases of more general formulae which will be derived in 
Section |]. 



14 



3.4 Using inversion relations to compute 
generating functions 

As in statistical mechanics, the inversion relation and symmetry, and some general assump- 
tions on analyticity of the generating function, are sometimes sufficient to determine the 
solution completely. In order to have an algorithm for computing a generating function term 
by term, it is necessary, but not sufficient, to have some property relating terms of different 
orders. For our purposes this property will always be x-y symmetry. Thus we restrict our 
attention to classes of graphs with x-y symmetry, i.e., Ferrers, staircase, directed convex, 
convex and three-choice polygons, and staircase polygons with a staircase hole. Moreover, 
we shall only consider the anisotropic perimeter generating function (without area). For 
the former four classes we will show that the inversion relation provides sufficient additional 
information to compute the generating function, whereas for the latter two it does not. 
The general form of the generating function is 

G(x,y) = H 1 (y)x + H 2 (y)x 2 + H 3 (y)x 3 + ■■■ (24) 

where the partial generating functions, H m (y) are rational functions 

HM = §§y (25) 

with D m (Q) = 1. The general form of the inversion relation is 

G(x, y) ± y a G(ex/y, l/y) = RHS (26) 

where a is an integer and RHS is zero or some simple function. It is equivalent to a self- 
reciprocity relation of the form 

H m (l/y)±e m y m - a H m (y) = RHS. 

Whether the inversion relation is sufficient to compute the generating function depends on 
the value of the exponent a and on the degree of the denominator, D m {y). Direct proof 
of the denominator form can often be obtained. For Ferrers graphs, it is easily shown that 
D m (y) — (1 — y) m - For staircase polygons, one finds D m (y) = (1 — y) 2m_1 . The same 
denominator form holds for directed convex and convex polygons also. For the three-choice 
polygons and staircase polygons with a staircase hole it can be shown that the denominators 
are 

,'(l-y) 2m ~ \l + y) 2m - 7 meven 
Dm{y) = \ (27) 

fl _ )2m-l( 1+ )2m-8 m odd . 
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Figure 4: A column-convex polygon 



We assume that in general we know the denominator form either empirically or by rigorous 
proof, and that D m (y) is of degree d m . 

Now we proceed inductively, following Baxter ||. If we have already computed the 
coefficient functions Hi(y), . . . ,H m -i(y) in the expansion (|2~4] ) and if x-y symmetry holds, 
we also know the coefficients of y,y 2 , . . . -,y m ~ l in the expansion of G(x,y). In particular, 
we can compute the coefficients of y, y 2 , . . . ,y m_1 in the numerator polynomial P m (y). In 
order to obtain the unknown coefficients of P m (y), we must be able to express them in 
terms of the known ones by means of the inversion relation. Writing P m {y) = Ylk a kU k i an d 
using D m (l/y) = ±y~ dm D m (y), the inversion relation fixes the value of the combinations of 
coefficients, ± ci£, with k + £ = a + d m — m. Hence the determination of all the coefficients 
at is possible if and only if the arithmetic condition holds: 



dm < 3m — a. 



(2? 



This condition is seen to hold for all the classes of convex polygons we have looked at, since 
d m < 2m — 1, but not for three-choice polygons or staircase polygons with a staircase hole, 
4m. 



since d m ~ 



4 Self-reciprocity via Temperley methodology 

We consider column-convex polygons as pairs of partially directed paths having the same 
endpoints, as indicated in Figure [|. 
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Let P be a column- convex polygon of width m. For < i < m, we denote by TV, 
(resp. Si) the number of north (resp. south) steps in the top path 7 at abscissa i. For 
< i < m, we denote by N_ t (resp. the number of north (resp. south) steps of the 
bottom path 7 at abscissa i. We choose the end points of the paths in such a way that 
jV = S = N m = S m = S = S m = 0. Note that 

m 

J2( N k + S k -S k -N k )=0. 



k=0 



We notice that all standard statistics are linear functions of the Ni^Si,^ and S^. For 
instance, the vertical perimeter of the polygon is 

rn 

2n = J2(N k + S k + S k + N k ) 

m 

= 2j2(N k + S k ). (29) 

k=0 

The height of the i th column of the polygon is, for 1 < i < m, 

i-i 

h t = J2(N k + S k -S k -N k ), 

k=0 

and the area of the polygon is 

m 

5>i - k)(N k + S k -S k -N k ). (30) 



m 

(I 

k=0 



Theorem 4.1. Let V be one of the following sets: Ferrers diagrams, stacks (drawn as in 
Figure \l(c} ), staircase polygons, bar-graphs, column-convex polygons. Let V m be the subset 



of V containing all polygons of width m. Let F m be the generating function for polygons in 
the set V m '- 

F m (y,z,y,z)= J] 

PeV m 

Then F m is a rational function, and it is self-reciprocal: 

F m (l/y, 1/z, l/y, 1/z) = C m F m (y,z,y,z), (31) 
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with 



' — m 



Vo 



| J y i for Ferrers graphs, 



i=l 



y2m—3 m—l 



Vo 



i=i 



m— 1 

- n ^ 

i=i 

:-ir 



y y m 



VoV r 



for staircase polygons (m > 2) 



for bar- graphs, 



for column- convex polygons. 



The proof of the theorem is based on the so-called Temperley approach for counting 
column-convex polygons [24], combined with the systematic use of formal power series [fF|. 
Here we provide only the proof for column-convex polygons, since the others are very similar. 

We commence by showing that the partial generating functions for column-convex poly- 
gons, V m (y,z,y,z), can be computed inductively. 

Proposition 4.2. Let V m (y,~z,y, z) be the generating function for column-convex polygons 
of width m. Let us denote it, for the sake of simplicity, V m (y ). Then the series V m (y ) 
can be defined inductively by: 



and 



v m +i(y r 



+ 



(i - y m Zm)( l - y m z m)v m (y m+1 ) 



-hn+l^JVmiZr 



+ 



i - y m+1 yJ(i - y^ +1 yj(zm - y n 



Proof. The basic idea is build a polygon of width m + 1 by adding a new column to 
a polygon of width m J7J. It is convenient to use Hadamard products to establish the 
functional equation. 
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Figure 5: Construction of column-convex polygon by Hadamard products. 



Let F(t) = J2 fht h and R(t) = E r ht h be two formal power series in t with coefficients 
in a ring A. We denote by F(t) R(t) the Hadamard product of F(t) and R(t), evaluated 
at t — 1: 

F(t)®R(t) = Y,hr h . 

In what follows, fh (resp. r^) will be the generating function for some column-convex polygons 
whose rightmost (resp. leftmost) column has height h, so that F(t)oR(t) will count polygons 
obtained by matching the rightmost column of a polygon of type F with the leftmost column 
of a polygon of type R. Also, R(t) will be a rational function of t. We shall use the following 
simple identity: 



F(t)Q 



1 — at 



F(a). 



The expression for Vi(y ) is obvious. We build a column-convex polygon of width m + 1 
as follows: we take a polygon of width m and match its rightmost column with the leftmost 
column of a column-convex polygon of width 2. 

This is illustrated by Figure |5|, which shows that 



V m {t)QR(t), 



(32) 



where 



R(t) 



ty , 



+ 



ty , 



1 — ty „ 1 — v ,«„ 1 — V 1 — ty „ 1 — y 1 — ty 



+ 



+ 



1 — ty , , 1 — £z m 1 — y , , 1 — ty 1 — tz m 1 — ty 
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is the generating function for column- convex polygons of width 2. In order to determine the 
coefficient of t h in R(t), we expand R(t) in partial fractions of t: 



R(t) 



2 



(l-^mK 1 -^) 1 

(1 — y ,y~r,)(l — y iZ m )(l — y~ l ,y )(l — y~ l _,z m ) 1 — ty 

(2m - y m+1 y m ^m) i 



(! - - y^ + i z ™){y m -z m ) i - 

+ 



l+l' 



(1 — y v^il — y 1 v )(z m — y ) 1 — ty 

Note that V^(0) = 0. We now combine eqn. ( [32] ) with the above expression for to 
obtain the announced expression for V m+ i(y m 1 ). a 



Proof of Theorem 4.1. Induction on m using the functional equation of Proposition 4.2 



shows that the partial generating functions for column-convex polygons satisfy: 

V m (l/y,l/z,l/y,l/z) = -^—V m {y,z,y,z). 

y y m 

■ 

We proceed similarly for the other families: the functional equation is obtained by setting 
some of the variables y~i,y.,Zi and z± to 0. Then, an inductive argument yields the self- 
reciprocity result. 

It would be tempting to write that the self-reciprocity of V m implies the self-reciprocity 
of, say, the generating function for staircase polygons, obtained by setting z { and Zi to 
in V m . But replacing a variable by in a self-reciprocal rational function might break the 
self-reciprocity: for instance, take P(yi, 2/2) = 1 + Vi + %i + 2j/ 2 + yxyi + y\yi- Then 

P(l/y h l/y 2 ) = J-P( Vl ,y 2 ), 
2/12/2 

but P(yi, 0) = 1 + yi + 2yf is not self-reciprocal. 

However, the following simple lemma gives a useful stability property of self-reciprocal 
rational functions. 

Lemma 4.3. Let F(yi, ■ ■ ■ ,y n ) be a self-reciprocal rational function. Let A be an m x n 
integer matrix. Let u = (ui, . . . ,u m ), and define u A to be the n-tuple whose i th coordinate 
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is Y\k uC k" ' ■ Then the series G{u) = F(u A ), if defined, is self-reciprocal in the variables Ui. 
More precisely, if F(l/y) = ±.y^F{y) then G(l/u) = ±u Af3 G(u) . 



From Theorem f4.1| and Lemma [4.3| we immediately deduce: 

Corollary 4.4. For any of the sets V listed in Theorem and any statistics on column- 
convex polygons that can be expressed as linear functions of the quantities N, S,N, S, the 
generating function for polygons in the set V m according to these statistics is a self -reciprocal 
rational function. 

This corollary allows us to complete the top part of Table |l|. Let us, for instance, derive 
the inversion relation satisfied by the tri-variate generating function G(x, y, q) for column- 
convex polygons, taking into account the usual parameters of interest: horizontal and vertical 
half-perimeters (variables x and y), and area (variable q). 

Eqns. fl29|) and (|30|) express the vertical perimeter and the area in terms of the quantities 
N, N, S and S. They imply that the (half) vertical perimeter and area generating function 
H m (y, q) for column-convex polygons of width m is 

H m {y,q) = V m (N,S,N,S) 
where yj. — z k — yq m ~ k and y k = z k = g~( m ~ fc ). Theorem JO] then gives 

H m (l/y,l/q) = l —H m (y,q), 

yq m 

which implies 

G(x,y,q)+yG(xq,l/y,l/q) = 0. (33) 

Note that in the first two self- reciprocity relations of Table 1, the exponent of q depends 
quadratically on the width. For this reason, they only yield an inversion relation for q = 1. 



5 Self-reciprocity via Stanley's general results 
5.1 Linear homogeneous diophantine systems 

Stanley has analyzed the situation where the objects to be counted correspond to integer 
solutions of a system of linear equations with integer coefficients (linear diophantine system) 
subject to a set of constraints. He has established certain conditions under which reciprocity 
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relations will hold between two combinatorics problems denned by the same linear diophan- 
tine system but by different sets of constraints, and also conditions under which the solution 
to a given problem will be self-reciprocal [^, [22|] . 

Consider the linear homogeneous diophantine system (LHD-system), 

= (34) 

in the unknowns ct = (cui, . . . , at s ) with <E» a matrix of integers having p rows and s columns 
and a p-tuple of zeros. The corank k of the system is defined to be s — rank($). For a 
linearly independent system, k = s — p. Let S be a set of integer solutions to eqn. ([34]). We 
define the generating function, S(y), as the formal power series 

s(y) = Y,y a (35) 

aes 

where y = (i/i, ... ,y s ) is a vector of fugacities associated with the unknowns in eqn. (3~1). 

In our applications, we find two types of constraints on the unknowns, ctj. Certain of 
the unknowns, ay, are required to be strictly positive while the rest are required to be non- 
negative. Conveniently, precisely these kinds of constraints have been treated by Stanley. 
Let the unknowns be ex. = (7, S) = 7 © S where 7 is an n-tuple and S is an (s — n)-tuple. 
Likewise let y = (it, v). In what follows, the notation, d > 0, means that all coordinates of 
S are positive. 

Proposition 5.1. Let E be the set of integer solutions, (7, 8), to a linear homogeneous 
diophantine system of corank k, such that 7 > and 6 > 0. Let E be the set of solutions 
to the same system with 7 > and 6 > 0. If the system has an integer solution, (7, S) 
such that 7 > and S < 0, then E(u, v) and E(u, v) are rational functions obeying the 
reciprocity relation 

E(u, v) = {-l) K E{l/u, 1/v). (36) 



Proof. This is Proposition 8.3 of ref. [20[ and the proof is given there. m 
Proposition |5.1| can be specialized to obtain a self-reciprocity condition, which will be 
our main tool in the derivations to follow. 

Corollary 5.2. A sufficient condition for the function E(u, v) to be self -reciprocal is that 
the linear homogeneous diophantine system has the solution (7, 5) = (1, —1). In this case 

u 1 

E(l/u,l/v) = (-ir—E(u,v). (37) 
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Figure 6: Staircase polygon of width three 



Proof. Since the solution (1,-1) satisfies the conditions of Proposition |5.1|, the reciprocity 
result fl36l ) holds. The result follows immediately from the shift (7, 8) — > (7 + 1, 5 — 1) which 
establishes a bijection between the sets E and E. u 

Since the conditions of the corollary are sufficient but not necessary, it is often possible 
to find a perfectly valid LHD-system describing a given self-reciprocal generating function, 
E(u,v), which does not admit the solution (1, —1). Hence we are faced with the problem 
of finding a suitable LHD-system which satisfies the corollary. A useful heuristic is to start 
with an LHD-system in many unknowns, and selectively eliminate those unknowns whose 
constraints are not independent of the constraints on the other unknowns. In all the cases 
we will consider, the resulting system will satisfy the conditions of Corollary |5.2| . We do 
not justify this heuristic here. In a paper subsequent to ref. pOf , Stanley |2l| develops a 
more comprehensive theory which overcomes these difficulties, and which additionally gives 
"correction" terms for systems in which self-reciprocity fails to hold. We have not yet 
explored the ramifications of this theory. 

Before applying the above result to staircase polygons with a staircase hole or to three- 
choice polygons, we use it to derive the reciprocity relation for ordinary staircase polygons 
of width three. This will serve to illustrate all the basic ingredients of the method. 



Example 5.3. Staircase polygons of width three can be characterized by the heights N±, 
N2, N_i, N^, Mq, Mi, M 2 and M 3 , as shown in Figure ||. Decomposing the polygon into 
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three columns, and imposing the condition that each column be as high on the left as it is 
on the right, we obtain the linear homogeneous diophantine system 

Mo - Mi - JVj = (38a) 
M 1 + N 1 -M 2 -N 2 = (38b) 
M 2 + N 2 -M 3 = 0. (38c) 

All heights must be nonnegative, but the self-avoidance condition additionally requires that 
the Mj be positive. The constraints M > and M 3 > are actually redundant, since they 
follow from eqns. ( |38a| , [38c| ) and the constraints on the remaining unknowns, namely 

N 1: N 2l N 1: N 2 >0 

M 1 ,M 2 >0. (39) 

Since the constraints on Mq and M3 play no role in the solution, we are free to eliminate these 
unknowns, and it turns out to be necessary to do so in order to apply Corollary |5.2|. We are 



left with the single equation ( |38b|) in the six independent unknowns 7 = (Ni , N 9 ,, N x , N ? ) 
and <5 = (Mi, M 2 ). Let us associate to the unknown N t (resp. N_ { , Mi) the fugacity y i (resp. 
V v Zi). 

Let E' be the set of solutions to eqn. ( |38b| ) subject to the the constraints 7 > and 
S > 0. Since 7 = 1, S = — 1 is a solution to eqn. ( |38b| ), Corollary [T^tells us that E'(y, y, z) 
is self-reciprocal, 

E'(l/y,l/y,l/z) = -^M^ E '(y,y,z). (40) 

ZiZ 2 

Equations ( |38a| , |38c| ) imply that to account for the dependent parameters M and M 3 , 
we make the substitutions z\ — > zqZi, y 1 — * zoy , z 2 — > z 2 z 3 and y 2 — > -23^2- Applying 
Lemma fO| , we obtain for the set E of nonnegative solutions to ( p8a|J38b| , |38c| ) such that 
M > 0: 

E(l/y, l/y, l/z) = - y -^^E(y, y, *). (41) 

Z\Z 2 

Notice that reintroducing the dependent unknowns has not changed the constant factor. 
This feature holds as well in the more complicated models we will look at. The result ( ffID 
may be verified by inspection of the explicit expression for the generating function 



E(y,y,z) 



z Z!Z 2 z 3 (l - y 1 z z 1 z 2 z 3 y 2 ) 



(1 - z y_^)(l - z 2i|/ 2 )(l - yi|/ 2 )(l - zqz x z 2 z- 3 ){\ - y x z 2 z 3 ){\ - y 2 z 3 ) 
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5.2 Applications 

We now apply the methods of Section ^TT] to staircase polygons with a staircase hole and to 
three-choice polygons. All the essential steps have already been seen in the derivation of the 
reciprocity result for staircase polygons of width three. They are 

1. Set up a linear homogeneous diophantine system by decomposing the polyomino into 
width one rectangles and imposing the condition that the left and right sides of each 
rectangle have equal height. 

2. Sort the unknowns into three classes, 7, d and r, according to whether they are 
constrained to be nonnegative, constrained to be positive or constrained by conditions 
on the other unknowns. 

3. Use Gaussian elimination to remove the unknowns in r. 

4. Verify that the resulting system is solved by setting all members of 7 equal to one and 
all members of 6 equal to minus one. Apply Corollary |5.2| to obtain the self-reciprocity 
result for the reduced system. 



5. Reintroduce the unknowns in the set r by means of Lemma [O 



We can define three widths for a staircase polygon with a staircase hole: the distance 
from the left edge of the figure to the left edge of the hole, k, the distance from the left edge 
of the figure to the right edge of the hole, £, and the width of the entire figure, m. Note 
that < k < £ < m. Recall that for staircase polygons the figures of width one were an 
exceptional case which did not obey the same reciprocity result as the general case. The 
staircase polygons with a hole of width one are also an exceptional case, which we must 
exclude. We thus impose the additional condition £ — k > 1. A figure with given k, £ and m 



is specified by the following dimensions, as shown in Figure 7TK 



1. heights N_j and Nj of the lower and upper perimeter segments of the polygon, 1 < j ' < 
m — 1, 

2. interior heights M~ to the left and right of, and within, the hole, < j < m, 



3. heights H_j and Hj of the lower and upper perimeter segments of the hole, k + 1 < j < 
£-1, 

4. interior heights M • and Mj below and above the hole, k < j < £. 
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(a) Staircase polygon with staircase hole, (b) Three-choice polygon, (k, £, m) = (2, 6, 5) 
(k,e,m) = (3,5,7) 

Figure 7: Labels for polyomino vertical heights 



Three-choice polygons can be regarded as staircase polygons with a hole which doesn't 
close. The width k has the same meaning as above, I denotes the ultimate horizontal extent 
of the branch of the figure above the hole, and m denotes the ultimate horizontal extent 
of the branch below the hole. Note that i > k and m > k. Again an exceptional case, 
m = k + 1, must be excluded. Hence we impose the restriction m > k + 1. The labeling 
of the vertical dimensions follows, with a few obvious modifications, the pattern of staircase 
polygons with a staircase hole and is shown in Figure |7(b)| . In particular, the heights Mj 
within the hole are defined only for j < mm(£, m). When £ = k the unknowns Mj and Hj 
do not appear. This special case is treated separately. 

As in the case of column-convex polygons, the standard statistics are linear in these 
heights. The (half-)vertical perimeter for staircase polygons with a staircase hole is given by 

m-l £-1 

n = M + J2^j + M k+ Yl **j ( 42 ) 

3=1 j=k+l 

and the area is given by 

fc-l £ £-1 m-l 

a = M + ^{Mj + Nj) + (Mi + M.j) + ^(Mj + Nj) + ^ (N^ + Mj) + M m . 

j=l j=k+l j=k j=£+l 

(43) 

In what follows, we associate to the unknowns Ni (resp. N_ { , Hi, H_ { , Mj, Mj, Mj the 
fugacities y { (resp. y., Wi, w i: Zi,z h z^. 

Proposition 5.4. Let Ek^ )m (y,y,w,w, z,~z, z) be the generating function for staircase poly- 
gons with a staircase hole where k, I and m are the widths defined above. Then if I — k > I, 
the generating function Ek^ m (y, y, w, w_, z, ~z, z) is self-reciprocal, 

E kAm (l/y, l/y, 1/w, l/w, l/z, l/z, l/z) = 

Zkzt 1^=1 {yjy) n5=l+i 



n7=i z iWj=k+iZjU 



— — E kAm (y,y,w,w,z,z,z). (44) 



j=l 1 lj=k+l -j 1 lj=k Z j 

Proof. The linear homogeneous diophantine system is the union of five sets of equations 
which we label Li-L 5 . The regions to the left and right of the hole give L\ and L 2 , the 
regions below and above the hole give L 3 and L 4 and the inside of the hole gives L 5 : 



M - Mi - Nj_ = 

Mj + Nj - M j+l - N j+1 = for 1 < j < k - 2 
M fc _! + iV fc _ 1 - M k -M k -M k -N k = 
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Mj, + M e + M t + N £ - M e+1 - Nj, +1 = 



Mj + Nj-Mj+t-N^-- 
M m _! + F m _i -M m = 







for I + 1 < j < m - 2 



M_ k -M k+l -N k+1 







Kj + Hj-Mj +1 -N j+1 



Mj + Nj - M j+ i - 
Mt-i + Nt^-Mt 



H 



M k - M k+1 - H 



■k+l 












for k + 1 < j < £ - 1 
for fc < j < t - 2 



(45) 



Mj + ^ - M j+1 - H j+1 = for k + 1 < j < £ - 2 
M,_x + - M, = 0. 

All heights of course are nonnegative. Self-avoidance imposes the additional constraint that 
the heights denoted Mj, Mj and M_j be positive. The set r, defined in step 2 above, contains 
six unknowns whose constraints are not independent which we eliminate as follows: M using 
the first equation of L 1; M m using the last equation of L 2 , M_ k using the first equation of L 3 , 
Mi using the last equation of L 4 , and M k and Mi using the first and last equations of L 5 . 
The resulting system is 



a 



Mj + Nj - M j+1 - N j+1 = for 1 < j < k - 2 
M fc _! + iVVi -M k - M k+1 - H k+1 - M k+1 - N k+1 - N k 

Mj + M e -! + He-! + M £ _ x + ivVi + N t - M e+1 - Ke +1 = 
Mj + Nj-Mj +1 -N j+1 



^ + H ] -M J+1 -N j+1 



M, 



Mj + Nj — ivi j + i — a j + i 



Mj + Hj-Mj +1 -H j+1 



-0 for £ + 1 < j < m - 2 

= for k + 1 < j < £ - 1 

= for k < j < £ - 2 

for A; + 1 < j < £ - 2. 



(46) 



The substitutions iVj, iV -, ifj, ^ 



1 and Mj, Mj, Mj 



-1 solve this new system of 
equations. One should note that when fc = lorm-l = 1 the system is somewhat modified, 
but one may check that the solution still holds. Therefore we may apply Corollary |5.2| 
to obtain a self-reciprocity condition on the generating function for the solutions of (J ■ L'j 



28 



subject to the positivity constraints on the heights. Making appropriate substitutions to 
restore the unknowns in set r, and using Lemma |4.3| we obtain eqn. 



We now treat three-choice polygons. 

Proposition 5.5. Let E kAm (y,y,~w,w,z,~z,z) be the generating function for three-choice 
polygons where k, £ and m are the widths defined above. Then if m — k > 1, the generating 
function E k> e )m (y, y, w, w, z, z, z) satisfies a self -reciprocity condition which, when £ = k, 
takes the form 

E KKm (l/y, l/y, l/w, l/w, 1/z, l/z, l/z) = 



ft— 1 — l \vn — 1 i fin — 1 

E k ,k,m(y,y,w,w,z,z,z), (47) 



n/c-1 — T—rm— 1 T—rm— 1 



Yij=i z j rij=fc+i £j 

and, when £ > k, takes the form 

E kAm (l/y, l/y, l/w, l/w, l/z, l/z, l/z) 



-E kAm (y,y,w,w,z,z,z). (Ah 



T-rmin(^m-l) TT^~1 -= TT m_1 y 
1 lj=l ^3 1 lj=k z j 1 lj=k+l -j 

Proof. It is simpler to treat the two cases £ = k and £ > k separately. The proofs follow 
very closely that of Proposition |5.4| . a 

As for column-convex polygons, the two propositions above may be extended to other 
statistics. 

Corollary 5.6. Let V be either of the sets staircase polygons with a staircase hole or three- 
choice polygons. Let Vk,£,m be the subset of figures in V with the widths k, £ and m defined 
as above. Then the generating function in V k ,i,m according to any statistics linear in the 
quantities N, N, H ; H, M ; M ; M ; is a self-reciprocal rational function (assuming £ — k > 1 
for staircase polygons with a staircase hole and m — k > 1 for three-choice polygons). 

The half-horizontal perimeter for either of the sets V is given by m + £ — k. Using this in 
combination with Corollary |5.6| , ( f4"2|) and fl43|) , we obtain the inversion relations specialized 
to horizontal and vertical perimeter, and area, which are listed in Table [3]. The exceptional 
cases (£ — k = 1 and m — k = 1 respectively) can be computed explicitly by the methods 
of 01. 
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6 Discussion 



Each of the methods we have discussed for obtaining reciprocity or inversion relations has 
its own particular uses. For example, the method of Stroganov is suitable for lattice models 
in statistical mechanics which are characterized by a family of commuting transfer matrices. 
The Temperley methodology is mainly applicable to families of polygons that are column- 
convex or nearly so. Stanley's method for obtaining reciprocity results apply to any problem 
defined by a system of linear homogeneous diophantine (LHD) equations, but the solutions 
to this system must be constrained by a system of simple inequalities of a certain form. 

It is probable that for many lattice models in statistical mechanics the low temperature 
expansion can be framed as an LHD-system. However, most are likely to require more 
general types of constraints than the simple inequalities of the directed polyomino problems 
we have considered. Likewise, the non-directed polygon problems that we have successfully 
treated using the Temperley methodology can be recast as LHD-systems with more complex 
constraints. How to handle such constraints is a worthy problem for future investigation. 

In recent work || this statistical mechanical language has been adapted for the enumer- 
ation of lattice paths, and may apply to polyomino problems as well. It is intriguing to 
speculate that the inversion relations found here may be connected with this approach. 

We have not searched for inversion relations for any polyomino problem in variables other 
than the natural variables for the problem. Yet the example of the Potts model demonstrates 
that such inversion relations may exist. It is also possible that symmetries in addition to 
the ones presented here can be found for some problems. It is our hope that such additional 
symmetries might lead to the solution of currently intractable problems. 

For the moment, we remark that the search for inversion and symmetry relations ap- 
pears to provide a new method to tackle certain combinatorial problems. The degree of 
applicability of this method is still unclear. 
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